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On a projection-corrected 

COMPONENT-BY-COMPONENT CONSTRUCTION 

Josef Dick*and Peter Kritzer^ 


Abstract 

The component-by-component construction is the standard method of finding 
good lattice rules or polynomial lattice rules for numerical integration. Several au¬ 
thors have reported that in numerical experiments the generating vector sometimes 
has repeated components. We study a variation of the classical component-by¬ 
component algorithm for the construction of lattice or polynomial lattice point sets 
where the components are forced to differ from each other. This avoids the problem 
of having projections where all quadrature points he on the main diagonal. Since 
the previous results on the worst-case error do not apply to this modified algorithm, 
we prove such an error bound here. We also discuss further restrictions on the choice 
of components in the component-by-component algorithm. 

Key words: Lattice point sets, polynomial lattice point sets, component-by-component 
algorithm. 
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Introduction 


Lattice point sets are integration node sets frequently used in quasi-Monte Carlo rules 
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for the approximation of s-dimensional integrals over the unit cube [0,1]®. For a modulus 
N {N a positive integer) and a generating vector g = {gi,... ,gs) G {1, 2,..., iV — 1}*, a 
(rank one) lattice point set is an integration node set of the form 



Here, for real numbers x > 0 we write {x} = x — [xj for the fractional part of x. For 
vectors x we apply {■} component-wise. A quasi-Monte Carlo rule using a lattice point 
set is called lattice rule. For further information on lattice rules we refer to OHSlEo]. 

*The research of J. Dick was supported under Australian Research Council’s Discovery Projects fund¬ 
ing scheme (project number DP150101770). 
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We consider a weighted Korobov space with general weights as stndied in [8l IT6] . 
Before we do so we need to introduce some notation. Let Z be the set of integers and 
let Z* = Z \ {0}. Furthermore, K denotes the set of positive integers. For a set £ 
we denote by \£\ the cardinality of £. For s G K we write [s] = {l,2,...,s}. For a 
vector X = (xi,..., G [0,1]^ and for u C [s] we write = {xj)j(zu G [0, and 
(a^u, 0) G [0,1]^ for the vector (|/i ,... ,ys) with yj = Xj if j G u and yj = 0 if j ^ u. 

The importance of the different components or groups of components of the functions 
from the Korobov space to be dehned is specihed by a sequence of positive weights 7 = 
( 7 n)MC[s], see [2T], where we may assume that 70 = 1. The smoothness is described by a 
parameter a > 1. The weighted Korobov space 'H{Ks^a,'y) is a reproducing kernel Hilbert 
space with kernel function of the form 


Ks,a,~j{x^y) 1 + 'y ^ 7ii y ^ 

0^«C[s] ^ 


|ii| 


exp(27ri/i„ • {x^ - yJ) 


It is well known in the theory of lattice rules that it is useful to restrict the range of a 
generating vector g of an Wpoint lattice point set to where 


Z^ = {ke{l,2,...,N-l}: gcd{k, N) = 1}. 


Furthermore, it is known (see, for example, |8]) that the squared worst-case error of a 
lattice rule generated by a generating vector g G Zf^ in the weighted Korobov space 
is given by 

e‘^ig)= ( 1 ) 

hGZ7{0} 
g h=0 (mod N) 

where for 0 7 ^ u C [s] and hu G Zi“' we have 

r( 7 , {h^,0)) = 7«n 

j&u 

It is known that the worst-case error in the Korobov space coincides with the worst- 
case error in the unanchored Sobolev space using the tent-transform [6] and is also related 
to the mean square worst-case error for randomly shifted lattice rules [5] . Hence the results 
here automatically also apply to those cases. 


The result 

The now standard method for hnding good generating vectors for numerical integration 
in Korobov spaces is the so-called component-by-component (CBC) construction (see [H 
m)- We can set the hrst component to 1 and then proceed inductively by choosing 
one new component at a time by minimizing the error criterion 92 , ■ ■ ■, Qd-i^ d) 

as a function of the last (not yet hxed) component g G Z^. Here, the components 
gl, 92 ,, g'^_i have been hxed in the previous steps. 

It has been observed that in running the CBC construction it may happen that com¬ 
ponents repeat themselves, i.e., there are i,j G {!,..., s} such that g* = g*. We quote 
from [T3] : 
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[.. .] However, it has been observed that the components start to repeat from 
some dimension onward for product-type weights, hence leading to a practical 
limit on the value of d [we remark that d has the role of s in This side 

effect of the CBC algorithm is yet to be fully understood. 

This problem may be due to numerical issues of the CBC algorithm, see [m p. 386], but 
this is currently not known. This paper also does not contribute to an understanding of 
this problem, instead we study a method to avoid its occurrence. Another quote is from 
P, where it is stated that: 

[.. . ] For large values of the worst-case error, the elements of the generating 
vector can repeat, leading to very bad projections in certain dimensions. 

To alleviate this problem, Gantner and Schwab [9] introduce a method they call pruning 
in the CBC algorithm. If gf ... ,gd-i have already been chosen by the CBC algorithm, 
then they choose the dth component from the set Zjsf \ {g[,... which forces the 

new component to differ from all the previous components. 

Following this idea, we study a modified CBC algorithm which excludes all values in 
a set Sd ^ Zn when choosing the dth component. Note that we allow the sets £d to 
depend on the values of ..., for d G {2,..., s}. (As just mentioned, [9] considered 
the special case £d = {gl, ■ ■ ■ The standard CBC algorithm can be obtained by 

setting £2 = ... = £s = i-e. no components are excluded in the CBC construction. We 
discuss other sets of exclusions in Section „ 

In the following, we write 0 for Euler’s totient function. 

Algorithm 1 Let A, s G N be given. 

(i) Set gl = 1, and choose £2 ^ Zjq. If no components are to be excluded in coordinate 
2, then set £2 = ^. 

(a) For d = 1,... ,s — 1, do the following: assume that gf...,g’^ have already been 
found, and choose £d+i ^ If no components are to be excluded in coordinate 
d + 1, then set £d+i = 0. Find as the minimizer g G Zj^ \ £d+i of 

In dimension d = 1 we do not use any exclusions since in our setting all components yield 
the same point set. 

One can still use the fast CBC method of [min] for this approach, with the additional 
step of checking whether a component is in the set of exclusions (as pointed out in the 
case of repeated components in [HI Section 4.2]). That is, in component d one needs 
to perform at most \£d\ checks for exclusions, hence one needs to perform at most an 
additional \£ 2 \ + ■ ■ ■ + \£s\ < (s — 1)(0(A) — 1) < sN checks. This does not increase the 
overall complexity of the fast CBC algorithm. 

Using exclusions in the component-by-component construction (for instance, by forcing 
new components to differ from the previous ones) implies that the theoretical results on 
the component-by-component construction as shown, for instance, in BEIE] do not 
apply anymore. Hence it has remained an open question as to what theoretical bounds 
one can get in this case. The following theorem provides an answer to this question. 
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For simplicity we assume in the following theorem that the weights are of the form 
7n = 7 j. However, it is clear that the result holds for any set of nonnegative numbers 

( 7 «)ucp]- We write C(tt) = Riemann zeta function. 


Theorem 1 Let N, s E N be given. Let 7 „ = Y\j^u 7i nonnegative real numbers 
7 j. Assume that g* = (^f^,..., ^f*) and sets of exclusions S 2 ,...,Ss ^ 2^ have been 
constructed by the algorithm. Then for all 1 < d < s we have 


r(gl,...,gS< 




l/A 


0 (iV) 


0 (iV) - \s, 


( 2 ) 


u<X[d\ jeu 

for any 1/a < \ <1, where the product over the empty set is defined as 1 . 


Proof. We prove the result by induction on d. For (i = l, letl/a<A<l. Then the 
result holds since 


e\9*i)=7{i} “ = iT2C(«)< (^2C(aA) 


l/A 


/igZ\{0} 


Assume that g* = (gl,... ,g'^) is chosen according to the algorithm and that ([ 2 ]) holds 
for d for any choice of 1/a < A < 1. 

From ([T]) it is easy to deduce that, for g E 

e\gl, ...,g:,g) = e\gl, ...,g:) + (3) 

where 

^iV,d+l,a,7(^) = ra{nr,h). (4) 

hgZ'i+ 1 , hd+i^O 
h-{g*,g)=0 (mod AT) 


From [H Eq. (5)] (setting (3j = 1) we obtain 

biN,d+l,a,~/i 9 ) 27d+lC(t^)A^ (1 + e {gi, ■ ■ ■ 1 gd)) T 'yd+lf^N,d+la,~/{g) 1 


where 

^N,d+l,a,~f{9) E E \hd+i\ "ro( 7 , h). (5) 

hd+ieZ h£7J^ 

N\hci+i h-g*=-hd+i9 (mod N) 

Hence we obtain from ([3|) 


e ( 9 ;,..., 93 ,g) = ( 1 + 27 j+iC(a)]V “)e ( 9 ;,... , 93 ) + 27 j+iC(a)]V “ +7J+iKjv,d+i, 0 , 7 ( 9 ). 

( 6 ) 

Let now A* E (1/a, 1] be chosen such that the right hand side of ([2]) for d + 1 is 
minimal for A*. In the following, we write 7 ^* for the weights 7^ = n,. u7j*- We apply 
Jensen’s inequality ® to obtain 


{e\gl,.. .,g:,g)r < {l+2^*7^l,aaX*)N-nie\gl, • • 


r-oA* 


*7aA* I oA*~A* 


r-oA* 
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Applying Jensen’s inequality to ([5]), we easily see that 


^ ^ (^Af,d+l,Q,7('^)) ^ ^ ^ ^Af,d+l,aA*,7^* (^) • ^N,d+l,a\* 




e^ZN 


ieZN 


In the following we use ideas similar to mi- We now use Markov’s inequality, which 
states that for a nonnegative random variable X and any real number c > 1 we have 
that P(X < cE(X)) > 1 — c~^. We use the normalized counting measure fi on as the 
probability measure. For c > 1 let 


Gc ■— {s' e Zjv : {KN,d+l,a,-f{g))^ < CKjv,d+l,aA*,7>* } 


— ^ S' ^ ■ {l^N,d+l,a,~f{9)) Z: ^ ^ {^N,d+l,a,'Y{^)) ^ • ^c- 


£eZjv 


Then 


m(Go) = ^ > m(Ac) = 


m) 




In other words, for any c > 1, there is a subset Gc C Z^ of size bigger than 0(iV)(l — c 
such that 

{f<-N,d+l,a,~f ig)r < c-Rn ,d+l,a\* ,'y^* 5 yg e Ge. 

By choosing c > 1 such that 




it follows that the set Gc \ 8d is not empty. This condition is satished for 

4>{N) 


c = 


0(iV) - \^d\ 


In particular, ii £d = ^ then c = 1. As is chosen by the algorithm such that the error 
e^{gl,... ,g^,g) is minimized, we obtain 


(e^gh • • ■,g*d,gki)r < (1 + 2^*7.";,C(«W)iV-“^*)(e^(^r,... ,^^))^* + + cfe 

< (1 + c2^*7.";iC(«A*)iV-“^*)(e2(^?*,.. .,g*d)r + c2^*7.";iC(«A*)iV-“^* + 

< (1 + c27,":,C(«A*)iV-“^*)(e2(^*,.. .,g*)r + c27,^;iC(«A*)iV-“^* + 


Using the induction assumption with A = A*, we obtain 

1 


ie^gl..., g*d)r < 


A A* 


0(iV) 


5^7r(2C(«A*))MJ] 


0(iV) 


HN)-\S,\- 


( 8 ) 


u£-[d] j£u 

Furthermore, from the proof of [U Lemma 5] we obtain 

< 9r(ry\*\(^ - lV-^^*\rh( r/T, X *'l 


Km+i,„a.,,»- S 2C(aA-)(l-]V-“"'),#.(]V)-' 7i'(2C(aA-))l 

07 ^uC[d] 


< 


2C(aA-)(l-]V-“"-).^(]V)-' 5 ^ 7 i-( 2 C(aA-))l“in 


MN) 


07 ^«C[d] 


jeu 


0(iV) - 


r(9) 
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Inserting the bounds in ([ 8 ]) and ([9]) into ([7]), and noting that N < (f){N) we obtain 

uQdl 

1 

0 

HN) 


jeu 


0(iV) - \s, 


+c2'y^_^_^({aX*) 


+c2'y^_^_^C{aX* 


0(iV) 

1 


0(iV) 


5 ; 7f(2C(aA*))l”IJ] 


07 ^nC[d] 


jeu 


0(iV) - |£, 


< 


cj>{N) 


7r(2CM*))HJ] 


0(iV) 


iiC[d+l] 


jGu 


0(iV) - |£, 


This implies the desired error bound for the special case of A*. However, since we chose 
A* such that the right hand side of (|2]) is minimal, the result also holds for arbitrary 
Ae(l/a, 1 ]. □ 

The following corollary considers the case where the relative size of the set of exclusions 
is uniformly bounded. 


Corollary 1 Let sequences of positive integers {Nk)km and {sk)km be given. Let g\ = 
{gfky--y9*k,k) be constructed by the algorithm using Nk,Sk and the set of exclusions 
S 2 ^\ ..., Ssp' for k Assume that there is a 0 < 6 < 1 such that 


I c(^) I 

md.X2<d<Sk \^y\ / r 

sup- 7 ' < 6. 

km <p[Alk) 


Then 


e\9l) < 





for any 1/a < X < 1 and any k eN. 


( 10 ) 


The corollary illustrates that as long as the relative size of the sets of exclusions is 
uniformly bounded ffTOj) . tractability results are not effected. In other words, if one, for 
instance, gets strong polynomial QMC tractability using the standard CBC algorithm, 
one also gets strong polynomial QMC tractability for the modihed CBC algorithm using 
uniformly bounded sets of exclusions. 


Remark 1 


(i) Similar results to the theorem and the corollary hold if one considers polynomial 
lattice rules (cf. instead of lattice rules (using an approach similar to instead 

ofW)- 

(a) Our method can also be applied to interlaced polynomial lattice rules flUUIfl- If the 
aim is to have different components in each coordinate, then due to the interlacing 
of consecutive coordinates in blocks of length d, it is enough to only force different 
first components of each block. This makes the additional construction cost and the 
increase in the error bound independent of the interlacing factor d. 
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Some particular choices for sets of exclusions 

Although the algorithm, the theorem and the corollary apply to arbitrary sets of exclu¬ 
sions, some particular choices are of general interest. We discuss some of them in the 
following. 

Repeated components 

If the aim is simply to avoid repeated components as observed in some numerical experi¬ 
ments, one can simply choose the sets of exclusions 

d = 2, ...,s. 

This makes sure that there are no two-dimensional projections of the integration lattice 
whose points all lie on the main diagonal {(x, x) : 0 < x < 1}. 

Avoiding diagonals 

To also exclude having two-dimensional projections where all the points of the integration 
lattice lie on an antidiagonal {(x, 1 — x) : 0 < x < 1}, one can additionally exclude the 
components N — gl, N — g^,N — g’^_i. This suggests to use the sets of exclusions 

Sd = {gl, N - gl, gl N-g;,..., g*_,, N - gU}, d = 2,...,s. (11) 

Note that this is only possible as long as 2(s — 1) < Even so, if, say 

for some small £ G N, then the factor (j){N)/— IT^I) = (j){N)/i becomes large, in 
which case the bound in the theorem becomes meaningless. So one still wants to impose a 
restriction of the form, say, maxi<d<s \£d\ < S(f){N) for some ‘reasonable choice’ (depending 
on the application) of 5 < 1. The last inequality implies that 2(s — 1) < 6(j){N). 

Avoiding diagonals in smaller dimension 

In some circumstances a condition of the form 2(s — 1) < 6(j){N) cannot be satisfied. For 
instance, when considering tractability questions one wants to study the dependence on 
the dimension as s tends to cxo [16]. Another case in which problems can arise is when N 
and s need to be increased simultaneously, as for instance in [T3] and |1]. In this case one 
can, for instance, choose Sd as in ([TTD for d = 1,..., s* for some fixed s* (independent of 
N and s) and set = 0 for d = s* -|- 1,..., s. As long as 2{s* — 1) < 5(p{N), the corollary 
applies since the relative size of the sets of exclusions is uniformly bounded and therefore 
strong polynomial QMC tractability results can be obtained. The particular choice of s* 
will depend on the problem under consideration. 

Reduced fast CBC 

Another instance where a particular type of sets of exclusions has been considered is 
[3|. In this case the aim was somewhat different, namely, to reduce the search space in 
coordinate d such that one obtains a speed-up of the fast CBC algorithm. In this case, 
instead of having to do additional computational work, the computational work actually 
decreases. See [3] for details. 
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